Orbital-ordering driven structural distortion in metallic SrCrOs. 



K.-W. Lcci^2 g^^^i Pickctti 
^Department of Physies, University of California, Davis, CA 95616, USA 
^Department of Display and Semiconductor Physics, 
Korea University, Jochiwon, Chungnam 339-700, Korea 
(Dated: October 1, 2009) 

In contrast to the previous reports that the divalent perovskite SrCrOa was believed to be cubic 
structure and nonmagnetic metal, recent measurements suggest coexistence of majority tetrago- 
nally distorted weak antiferromagnetic phase and minority nonmagnetic cubic phase. Within the 
local (spin) density approximation (L(S)DA) our calculations confirm that a slightly tetragonally 
distorted phase indeed is energetically favored. Using the correlated band theory method (LDA+ 
Hubbard U) as seems to be justified by the unusual behavior observed in SrCrOa, above the crit- 
ical value Uc=4: eV only the distorted phase undergoes an orbital-ordering transition, resulting in 
^2g ~^ dl,y{dxzdyzY corresponding to the filling of the d^y orbital but leaving the other two de- 
generate. The Fermi surfaces of the cubic phase are simple with nesting features, although the 
nesting wavevectors do not correlate with known data. This is not uncommon in perovskites; the 
strongly directional d — d bonding often leads to box-like Fermi surfaces, and either the nesting is 
not strong enough, or the matrix elements are not large enough, to promote instabilities. Fixed spin 
moment calculations indicate the cubic structure is just beyond a ferromagnetic Stoner instability 
{IN{0) ~1.1) in L(S)DA, and that the energy is unusually weakly dependent on the moment out 
to l.S/is/Cr (varying only by 11 meV/Cr), refiecting low energy long-wavelength magnetic fiuctu- 
ations. We observe that this system shows strong magneto-phonon coupling (change in Cr local 
moment is ~7.3 /^s/A) for breathing phonon modes. 



PACS numbers: 71.20.Bc, 71.30.-|-h, 75.50.Ee 



I. INTRODUCTION 

Forty years ago, a few divalent chromate per- 
ovskites ^CrOa (^=Pb, Sr, Ca), formally possess- 
ing the Cr^"*" ion, were synthesized at high temper- 
ature ~1300 K and under high pressure 6-10 GPa 
by a few groups. In spite of their atypi- 
cal and controversial properties, these systems have 
been little studied, probably due to difficulty of syn- 
thesis. More recently a few groups have begun to 
revisit the CaCrOa and SrCrOa compounds. @, 0, 

H, Whether these systems are metallic, strongly 
correlated, and spin ordered is still controversial. [6|, 

UMMM 

Roth and DeVries reported an ordered moment of 

I. 9 and Curie-Weiss moment of 2.83 hb in the 
isovalent compound PbCrOs, consistent with S = 1 
Cr^+ (rf^).[H Chamber land and Moeller synthesized 
a single crystal, which was semiconducting with 0.27 
eV activation energy. There was an anomaly at 
Ti=:240 K, and an upturn at T2=160 K, in suscep- 
tibility. The latter was thought to imply a G-typc 
antiferromagnetic (AFM) ordering corresponding to 
antiparallel spin ordering between all nearest neigh- 
bor Cr'*+ ions. Additionally, at T3=100 K, the log- 
arithmic resistivity shows a kink, implying another 
transition. The samples of both groups had a cu- 
bic structure with lattice constant a w 4.00 A. Lo- 
cal (spin) density approximation (L(S)DA) calcula- 
tions obtained a magnetic moment of 1.4 /xs, three- 



quarters of the experimental value, but no band- 
gap. Ill This difference points to interaction effects 
beyond those described by L(S)DA. Chambcrland 
also synthesized a cubic SrCrOa with a=3.818 A, 
and concluded it to be a nonmag netic (NM) metal. Q 

Goodenough et al. synthesized polycrystalline 
CaCrOa, which is orthorhombic (a=5.287 A, 
5=5.316 A, and c=7.486 A) and non-conducting (al- 
though probably due to polycrystallinity) . Q Wei- 
her, Chamberland, and Gillson obtained a metallic 
single crystal sample. Q The susceptibility measure- 
ments showed two anomalies, a kink at 325 K and 
an upturn at 90 K. At the latter, which is recently 
identified as the Neel temperature TAr,0, a kink 
in the resistivity data and decrease in volume by 
2 % from the powder diffraction studies were ob- 
served. The Curie- Weiss moment is high spin 3.7 
fiB, recently confirmed by Zhou et aZ.Q The reason 
of large difference from spin-only value of 2.8 ^b for 
S—l system is unresolved, however it was found that 
the susceptibility in SrCrOa did not follow a Curie- 
Weiss behavior so no local moment value could be 
identified. 

In more recent studies, Zhou et al. observed 
a smooth decrease in thermal conductivity of 
CaCrOa and SrCrOa compounds as temperature is 
lowered, in their interpretation characteristic of nei- 
ther an insulator nor a metal. Q They interpreted 
this unusual behavior as due to some unusual Cr- 
O bonding instability, supported by an increase in 



compressibility observed around 4 GPa. In contrast 
to Chambcrland's initial suggestion. Q Zhou et al. 
concluded that SrCrOs is NM and insulating. Ko- 
marck et al. reported antiferromagnetism with or- 
dering wavevector Qm = (|,|,0) (in units of ^) 
and a saturation magnetic moment of 1.2 /i^, using 
SQUID susceptibility and neutron diffraction. Q At 
Tjv, in contrast to the preliminary observations, § 
no change in volume was apparent. They suggested 
CaCrOs is itinerant, but close to being localized, 
implying importance of correlation effects. [9l.lToj Ad- 
ditionally, no evidence of orbital ordering within the 
t2g shell was observed. 

Attfield and coworkers have concentrated on 
SrCrOa, using neutron diffraction and synchrotron 
powder x-ray diffraction studies. [7i,i8i] Below Tm ~40 
K, most of their sample underwent a structure tran- 
sition from a NM cubic phase to a AFM tetragonal 
phase, but the two phases coexist even at low T. 
In the tetragonal phase with only slightly inequiv- 
alent lattice parameters, at T/v there are no visi- 
ble changes in volume or in averaged Cr-0 distance, 
and additionally no kink in resistivity. The tem- 
perature dependent neutron diffraction data imply 
orbital reoccupation {dxyd^zdyzY -^dl.y{dxzdyzY 
the tetragonal phase at or near Tjq. Both phases 
are metallic, though showing high resistivity due to 
grain boundary scattering. 

In this paper, we will focus on SrCrOs, which 
has been little studied theoretically to date. In Sec. 
Ill, in particular, we will address electronic struc- 
tures of both NM cubic phase and AFM distorted 
phase, and show that the inclusion of correlation ef- 
fects leads to an orbital-ordering driven distortion 
while remaining metallic. In Sec. IV, we discuss the 
oxygen breathing vibration that shows quite strong 
magneto-phonon coupling. 



II. STRUCTURE AND CALCULATION 

Attfield and coworkers suggested a small struc- 
ture distortion involving relative displacement of Sr 
and O ions, leading to ^/2a x ^/2a x 2a quadru- 
pled supercell (space group: Imma, No. 74). Q 
The corresponding AFM order we consider, which is 
(i, i, i) order in terms of the original perovskite cell, 
is a G-type AFM. Through more precise measure- 
ments, instead of the quadrupled supercell they more 
recently concluded a tetragonally distorted struc- 
ture with c/a saO.992 but negligible change in vol- 
ume (less than 1%), which shows antifcrromagnctic 
ordering at wavevector Qa/=(|,^,0) (space group: 
PA/mmm, No. 123). Q This AFM order is conven- 
tionally called C-type AFM. Although this tetrago- 
nal phase is dominant below Tjv, a NM cubic phase 
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FIG. 1: (Color online) Enlarged band structure of NM 
SrCrOa in the region having mostly Cr d character. A 
flat band along theV-X -M hne hes at 0.5 eV. The O 
2p states lie on the regime of -7.5 to -1.8 eV (not shown 
here). These symmetry points follow a simple cubic no- 
tation (see Fig. [3]). The R point is a zone boundary along 
(111) direction. The horizontal dashed line indicates the 
Fermi energy Ep. 



also coexists, consistent with results of our calcula- 
tions (see below). 

We investigated both distorted structures as 
well as the cubic phase, using L(S)DA and the 
LDA-^Hubbard U (LDA-I-U) method. [H, [H For 
the quadrupled supercell, only the planar O position 
is displaced along the (001) direction, since negligi- 
ble displacements in Sr and apical O were initially 
suggested. We used the recent experimental lattice 
constant a=3.811 A for the cubic structure and for 
the quadrupled supercell, and a=3.822 and c=3.792 
A for the tetragonal structure. Q Our optimized lat- 
tice constants in the cubic phase within L(S)DA are 
3.748 A for NM and 3.76 A for AFM, reflecting 
the usual small increase in volume associated with 
magnetism and the overbinding that is common in 
L(S)DA. 

Wc have used the full-potential local orbital code 
FFLO for our study. In FPL0-5,[ll basis orbitals 
were chosen as Cr (3s3p)4s4p3d, Sr (4s4p)5s5p4d, 
and O 2s2pid. (The orbitals in parentheses indi- 
cate semicore orbitals.) The Brillouin zone was sam- 
pled with a regular mesh containing 726 irreducible 
k points, since a fine mesh is required for sampling 
the Fermi surface. 



2 







Total ! 




Cr 1 j 1 




-—- o i j i 




PMEy I i 




■ FMEp ! '\j 




----- FME^ j / 


\ 


i /' 

1 \ J ,1 


\\ 
\\ 




\ /^'^^ ! 


~ — — 



I 1 ■ ^ 1- - 7 - . I ,_, 

-2-10 1 

E-E^ (eV) 



FIG. 2; (Color online) Total and atom-projected densi- 
ties of states in NM SrCrOs. A sharp peak lies at 0.5 
eV, while around Ef there is only a smoothly monotonia 
behavior. The DOS at Ef N(0) is 1.53 states per eV per 
spin, a quarter of which is contributed by O ions. 



III. ELECTRONIC STRUCTURE 
A. NM Cubic Phase within LDA 

In this subsection we address the electronic struc- 
ture of the NM cubic phase (observed above 40 K) 
using LDA. The band structure around the Fermi 
level Fip (the Cr d regime) is shown in Fig. [l] and 
the corresponding densities of states (DOSs) are dis- 
played in Fig. [21 The one-third filled Cr t2g manifold 
with width of 2 eV lies between -1.5 eV and 0.5 eV 
(we take Ej? as the zero of energy). The unfilled Cr 
eg manifold touches the t2g manifold at the X point 
at 0.5 eV and extends to 4.5 eV, leading to the t2g- 
Cg (midpoint) splitting of roughly 2.5 eV. Flat bands 
along the T — X ~ M lines result in a sharp peak in 
the DOS at 0.5 eV, which otherwise does not have 
any distinguishing structures near Ej;^. 

The Fermi surfaces shown in Fig. [3]display nesting 
features, indicative of large susceptibilities at related 
wavevectors and suggesting the possibility of either 
magnetic or charge instabilities. The intersecting 
pipe-like surface has six circular faces with a radius 
of 0.32(-^), so contains ~0.6 electrons per spin. Two 
cube-like surfaces with sides of length ~0.61(^) are 
very similar in size, and touch along the T — X line. 



B. Structure Relaxation 




FIG. 3: (Color online) Fermi surfaces, which contain 
electrons of NM SrCrOa. These surfaces show nesting 
features commonly arising in perovskite structure tran- 
sition metal oxides. Both the second and the third sur- 
faces are cube-like with rounded edges and corners; each 
face of the second surface is nearly circular. The Fermi 
velocity of 2x10^ cm/sec, which is a typical value in a 
metal, is nearly uniform through the surfaces. 



Our L(S)DA calculations show this tetragonally 
distorted structure is energetically favored over the 
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FIG. 4: (Color online) Changes in energy SrCrOs with 
respect to lattice parameter ratio (c/a) in C-AFM phase, 
with various choices of lattice parameter a. This a pa- 
rameter in the range of 3.70-3.82 A is approximately the 
range between our optimized value in the cubic structure 
and the experimentally observed value. The arrow de- 
notes the change in position of the minima, pointing out 
the stronger structural distortion for larger volume. 

cubic structure. As expected from the small change 
in structure, however, the difference in energy is 
small, no more than 1 meV per Cr. The fact that 
these two structures are nearly degenerate is consis- 
tent with experimental observations. The quadru- 
pled supercell, on the other hand, has a slightly 
higher energy than the cubic phase. In this sec- 
tion, we will focus on only the competing cubic and 
tetragonal phases. 

To investigate sensitivity of this structure distor- 
tion to volume, we calculated the energy vs. c/a re- 
lation in the range of a=3.70-3.82 A. The volume is 
kept fixed while the c/a ratio of lattice parameters is 
varied, since the experiment shows negligible change 
in volume between the two observed phases. Q The 
result is given in Fig. |4l In very close agreement 
with experimental observations, a minimum occurs 
at c/a=0.99 for a=3.82 A, which is the experimen- 
tally observed lattice parameter. However, with de- 
creasing volume this distortion is gradually relieved, 
and finally the cubic structure is favored energeti- 
cally somewhat below a ~3.70 A, 0.06 A smaller 
value than our optimized parameter. 



C. Energetics 

As discussed in the Introduction, the recent obser- 
vations indicate the coexistence of NM cubic phase 
and AFM tetragonally distorted phase. Q Within 
L(S)DA, AFM order is more favored energetically 



FIG. 5: (Color online) Fixed spin moment calculations 
using both one Cr and two Cr cells. The latter cell allows 
both NM and AFM states at moment M=0. For the 
doubled cell, the plot shows richer behavior (see text). 
Above IjiB the states are the same simple FM aligned 
state. 



than NM in both crystal structures. In the tetrag- 
onal phase, AFM order (M=1.55/x_b) is favored en- 
ergetically over FM (M=1.16/iB) by 150 meV/Cr. 
The ferromagnetic state is favored only 11 meV/Cr 
over the nonmagnetic state. With Stoncr / = 0.6 
eV (see below), the energy gain in a simple Stoner 
picture /Af^/4 would lead to much larger value of 
^0.2 eV/Cr. Note that the moments are roughly 
consistent with S=l ions in the presence of strong 
p^d hybridization that is evident from the projected 
DOS in Fig. [21 The strong reduction in moments 
from the ionic value also reflects considerable hy- 
bridization, so a fixed-local-moment (i.e. Heisenberg 
picture) has little use here. The behavior of ferro- 
magnetic order will be investigated more explicitly 
below by fixed spin moment calculations, [l^ The en- 
ergy differences between magnetic phases are nearly 
independent of the small distortion that we study 
here. 



D. Fixed Spin Moment (FSM) Studies 

Initially we used a single Cr cell for our calcu- 
lations, supplementing this with doubled cells (see 
below). Consistent with our L(S)DA calculations, 
E(M) has a minimum at M ~1.3 /is and a related 
small gain of 11 mcV in energy, as can be seen in 
curve in Fig. \5\ The very flat E(M) behavior for M 
up to l.bfiB indicates that magnetism in SrCrOa is 
very peculiar. The gain in exchange energy IM"^ /A 
is almost exactly compensated by a cost in other en- 
ergy contributions across this range. The energy vs. 
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moment curve is fit at small M to the expression 
e — eo = aM^ + /3M^ to evaluate the Stoner (ex- 
change) constant / from these FSM calculations. |15| 
The resulting value of a = —17 meV//x^ (and /3 = 6 
meV//i^) provides the enhanced (observed) suscep- 
tibility given by 



Xo 



1 - N{Q)I 



(1) 



where the bare susceptibility is xo = 2^^-/V(0). The 
Stoner enhancement factor S = [2axo\~^ is about 
-8. Thus the Stoner / = 0.6 cV, and IN{Q) = 
1.1 with A''(0) = 1.87 states per eV per spin from 
our FM calculations, predicting the system is beyond 
the Stoner magnetic instability in the cubic phase 
(within L(S)DA). 

To generalize the study, we used the two-Cr super- 
cell and started from C-AFM order (G-AFM order 
showed similar change) which also has total moment 
M=0. The magnetic field that is applied in the fixed 
spin moment method provides an evolution of the 
C-AFM state into a FM state, perhaps through an 
intermediate ferrimagnctic phase. A small system- 
atic energy difference between the single and dou- 
bled cell energies that we compare has been accom- 
modated by aligning the M=0 energies. We find the 
energy vs. M curve to be comprised of three sep- 
arate regimes: AFM at M = where the energy 
is ^^150 meV/Cr lower than for NM; ferrimagnctic 
for < M < 1 where the two Cr moments differ; 
then FM for M > 1 where both calculations de- 
scribe the same simple FM phase. While one might 
expect strong magneto-elastic coupling in this sys- 
tem, we find that our FSM results are insensitive to 
c/a ratio in the range given in Fig. |3) 



E. L(S)DA Electronic Structure of the AFM 
Tetragonally Distorted Phase 



Although Attficld and coworkers observed coexis- 
tence of the nonmagnetic cubic phase and the an- 
tiferromagnetic tetragonal phase, our L(S)DA cal- 
culations show energetically favored AFM in both 
phases, as already addressed. The C-AFM band 
structure is shown in Fig. [S] The Cr local moment 
is 1.55 yUs, with negligible dependence on this small 
tetragonal distortion. The largest effect of the struc- 
ture distortion on the band structure occurs in the 
maxima at the M and R points, downshifting in en- 
ergy at most 25 meV. The topmost band crossing 
Ep has the dxy character. 




FIG. 6: (Color online) Enlarged band structure of C- 
AFM in the tetragonally distorted structure, plotted in 
the basal plane. That of the cubic phase with the same 
volume is similar, so is not shown. The antiferromag- 
netism introduces a gap at the X point in the range of 
-0.5-0.5 eV, and the t2g bands become disconnected from 
the higher eg bands due to the reduced bandwidths. In 
unit of vr/a, the X point is (1/2,1/2,0). 
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FIG. 7: (Color online) [/-dependent orbital-projected 
densities of states of Cr t2g states for (a)-(b) the cubic 
and (c)-(d) the tetragonal phase, at U=0 and 4 eV in 
C-AFM order. At U—4 eV, an orbit-ordering transition 
in dxy orbital occurs in the tetragonal phase. 



IV. INCLUSION OF CORRELATION 
EFFECTS 



L(S)DA predicts the C-AFM phase to be consider- 
ably lower in energy than the NM phase, in disagree- 
ment with experimental data that suggests the two 
phases are nearly degenerate. Also L(S)DA cannot 
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produce the partially orbital order experimentally 
suggested. We address this discrepancy by includ- 
ing correlation within the LDA+U approach. On- 
site Coulomb repulsion U was applied on the Cr ions 
with AFM order. U was varied in the range of 0- 
8 eV, but the Hund's exchange integral J=l eV is 
fixed since the results in the physical range of U are 
expected to be insensitive to J. 

At U=0=J (i.e., LSDA level), two electrons are 
evenly distributed in the three bands of the major- 
ity t2g manifold. In the C-AFM (tetragonal) phase, 
increasing U changes occupancies, with d^^ and dy^ 
weight transferring into dxy At C/c=4 eV, the ma- 
jority dxy band is fully occupied, while the majority 
dxz and dyz bands share equally the other electron, 
as shown in the orbital-projected DOS given in Fig. 
[71 This dl.y{dxzdyzY state is consistent with Attfield 
and coworkers' conclusion from their neutron diffrac- 
tion measurements.^] This partial orbitally ordered 
arrangement remains for larger values oi U . A po- 
tentially Mott-insulating dl.^,dy^,d!^y state is avail- 
able but did not arise in the calculations. In the 
range of U studied here, this system remains metal- 
lic. 

An orbital ordering transition in a multi- 
band system can result in significantly different 
bandwidths.[T6| However, in this case, the differ- 
ence in the occupied bandwidths is small, about 100 
meV. The important change is that the center of 
the dxy band lies roughly 1 eV lower than those of 
the partially occupied dxz,dyz bands. It is this dif- 
ference in the band center (i.e. the on-site energy) 
that drives this transition, as happened in Naa;Co02 
p7j or V2 03.[l3| Keeping the structure (and sym- 
metry) cubic inhibits such an orbital-ordering tran- 
sition, implying close interplay between structural 
distortion (even though tiny) and orbital ordering 
in LDA+U calculations in this system as suggested 
in CaCrOg.H 



V. O TILTING AND BREATHING 
PHONON MODES 

To investigate another possibility of structure dis- 
tortion in this perovskite, we used G-AFM order 
in the quadrupled supercell, which allows O tilting 
and breathing phonon modes in frozen phonon cal- 
culations within L(S)DA. However, our calculations 
show these distortions are unfavored energetically, 
resulting in stable phonon modes. These data of 
change in energy vs. O displacement are fit well with 
a simple harmonic function. First, for the O tilting 
vibration the phonon energy is 38 meV, typical for 
metallic oxides. For the breathing mode the energy 
is 89 meV, corresponding to an rms displacement 
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FIG. 8: (Color online) Change in moments due to the 
O breathing vibration when AFM allows. The Crl, 
Cr2, and total moments are changed 7.1, —7.5, and 15.3 
^fl/A in magnitude. Note that no displacement leads 
to AFM. 



of the oxygen ions by 0.05 A. (Allowing magnetic 
ordering, both frequencies reduce by ^ 5%.) 

Our calculations show strong magneto-phonon 
coupling for the breathing mode. When AFM or- 
dering is included, the Cr local moment of Mcr= 1.5 
[iB is modulated by about ±7.3 ^b/^^ as shown 
in Fig. [SI These changes are quite large, even 
larger than the change in Fe moment, 6.8 /is/A, 
in LaFeAsO when As ions are displaced, which 
is widely discussed as unusually strong magneto- 
phonon coupling. (l9| At the rms displacement, the 
difference between Cr moments becomes 0.73 /^s, il- 
lustrating just how large the modulation of the mo- 
ment by the breathing mode is. The Cr charge "dis- 
proportionation" due to O breathing also shows a 
value of ±0.13e {i.e. a 0.26e charge difference) at 
the rms displacement, corresponding to a shift in 
charge of ±2.6e/A. 



VI. SUMMARY 

We have presented and analyzed the electronic 
structure, magnetic ordering, and the impact of 
strong correlation effects in the perovskite material 
SrCrOa, which is reported in a nonmagnetic cubic 
phase coexisting with an antiferromagnetic phase 
concurrent with a small distortion in structure below 
40 K. Consistent with these observations, L(S)DA 
predicts the slightly distorted tetragonal structure, 
but overestimates the polarization energy. With 
L(S)DA the cubic phase is magnetically unstable, 
with a Stoner product IN(0) «1.1. The Fermi sur- 
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face shows nesting features, but they do not cor- 
relate with known data. Although it is possible for 
certain crystal symmetries for orbital orderin g to oc- 
cur withouth changing the crystal structure, [lOj we 
have not pursued such possibilities. 

Including correlated effects within the LDA-I-U 
approach, the distorted phase undergoes an orbital 
ordering transition at a critical interaction strength 
Uc=4: eV, leading to ^ dl.y{dxzdyzY orbital or- 
dering and structural transition to tetragonal struc- 
ture. The structural symmetry lowering is crucial; 
the cubic phase remains a simple metal even for 
higher U . We have also demonstrated that the O 
breathing modes show strong magneto-phonon cou- 
pling. 

The magnetic behavior in SrCrOa remains un- 
clear. The experimental data show weak magnetic 
behavior, and the susceptibility does not follow 
Curie- Weiss behavior. @ In contrast to these obser- 
vations, our various calculations always result in a 



full moment corresponding to S=\ configuration (re- 
duced by hybridization) as expected for a ion. 
Without any peak in the DOS, the temperature vari- 
ation cannot be modeled with temperature broaden- 
ing as can be done, for example, in TiBe9.|i20i] 
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